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Chapter 1 
Introduction 


This progress report supports the development of a bistatic radar cross 
section (bistatic RCS) measurement range by PMTC Milicon-P986 under 
NASA Lewis Research Center contract #NAG3-1125. The bistatic scat- 
tering measurement scenario under consideration is shown in Figure 1. As 
shown in this figure, a compact range consisting of a parabolic reflector and 
a spherical wave feed will be used to generate a plane wave that illuminates 
the target. The incident plane wave signal is bistatically scattered from 
the target under test. The bistatic RCS measurement system must then 
determine the far-field scattered signal in the bistatic direction. 

Many methods exist to estimate the far-field scattered signal without the 
distances necessary to measure the actual far-fields. One frequently used 
method is the reception of a plane wave by a compact range reflector. Other 
methods include using a one-dimensional vertical parabolic panel, such as 
a bent rod, scanned in a radial manner around the target. This measure- 
ment method requires a one dimensional near-field to far-field transforma- 
tion algorithm. Another alternative is to scan the scattered fields along a 
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Figure 1: Bistatic Scattering Arrangement 
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spherical, cylindrical or planar surface with a probe antenna and transform 
the measurements to the far-field. Each of these surface scans require a 
two dimensional transformation algorithm, but planar scanning does not 
require curved sampling paths and it uses a well known near-field to far- 
field transformation algorithm, the heart of which is a two-dimensional 
Fourier transform [2]. A fourth alternative is synthesis of a plane wave by 
an antenna array, recreates the aperture distribution of a compact range re- 
flector. In considering each of these alternatives it is necessary to determine 
the effects of finite scan area and sample spacing. 

Note the similarity of these last two alternatives to the classical planar 
near-field technique for antenna pattern measurement. Upon illumination 
of the scatter by the plane wave, the scatterer becomes a radiator. It 
is behaving like an antenna. This means that much of the planar near- 
field technology developed for antenna pattern measurements is applicable 
to this research. For the measurement of antenna patterns, knowledge of 
the gain of the antenna is required such that the pattern function can 
be normalized to the total radiated energy. Thus if 85% of the energy 
is received by the scan surface, the gain can be well estimated. For the 
measurement of RCS however, only the plane wave incident at the given 
bistatic (or monostatic) angle is required. This is accomplished by any 
planar aperture large enough to reduce edge effect errors. Therefore, a RCS 
measurement system need not receive 80 to 85% of the radiated energy from 
the target. A RCS measurement system is calibrated by the theoretical RCS 
predictions of known test targets. 

In this progress report, planar near-field scanning will be considered. 
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The scan plane with reference coordinates is shown in Figure 2. Reciprocity 
indicates that if an antenna array can simulate a compact range antenna 
pattern, then the antenna array should also receive signals similarly to a 
compact range reflector. 
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Figure 2: Scan Plane Reference Coordinates 
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Chapter 2 

Scanner Design Issues 

2.1 Introduction 

Scan plane size (array size), probe measurement spacing (array element 
spacing), position accuracy, probe characteristics (array element character- 
istics), and the distance from the scanner (or array) to the test target all 
effect the scanner system and an antenna array system similarly during the 
design stage. 

2.2 Scan Plane Size 

In classical near-field antenna pattern measurements, near-field to far-field 
transformations, such as Kern’s Plane Wave Spectrum matrix approach [10], 
direct determination of the source distribution at the target, or the Lorentz 
reciprocity formulation developed by Paris, Leach and Joy [7] can be used 
to estimate the far-field scattered signal from the measured scan plane data. 
To quantize the magnitude of the error involved in these transformations 
over a finite scan plane, Yaghjian [8] suggests performing the transformation 
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algorithm on the area outside the scan plane. Since the fields outside the 
scan plane are not known, the error is estimated as the “a priori” knowledge 
of the energy not received by the scan plane. This criterion is frequently 
used to require retaining 80 to 85% of the radiated energy within the scan 
surface and limit test targets to ones that radiate directively. 

In this work, the array size and spatial windowing designs will be cho- 
sen to mimic the behavior of a compact range reflector. The array size 
requirement is thus the same as for a compact range reflector: the plane 
wave region must cover the test target entirely. For identical target (quiet) 
zones, the necessary scan plane size is expected to be comparable with a 
compact reflector surface (including rolled edges). 

The size of the plane wave region in the target zone is proportional to the 
size of the scan plane. A useful plane wave region is often characterized by a 
1 dB amplitude ripple and 10 degree phase ripple. In the literature, Johnson 
et.al. [1] quotes Woonton as having a minimum scan-plane to target-zone 
size ratio of 2.0, Chapman 2.3, and Georgia Tech 1.9. Joy and Rowland [3] 
used a minimum scan-plane to target-zone ratio of 1.7 for a highly directive 
parabolic dish target. Thus, as a general rule of thumb, the scan-plane must 
be twice the size of the largest target to be measured. 

2.3 Probe Measurement Spacing 

The planar scanner receives the incident scattered fields at discrete sample 
points. To satisfy the Nyquist sampling criterion, samples should be taken 
at intervals which are equal to or less than one half the wavelength. Cown 
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and Ryan [9] however, note that “thinned” near-field data sets can be used 
to find accurate RCS values over reduced angular sectors. The grating lobes 
that occur due to the increased sample spacing are outside the region of 
interest and are simply ignored. If the incident planar-wavefront angles are 
restricted to a small range of angles about the normal of the scan surface, 
the sample spacing can be increased beyond the half-wavelength criteria. 

The maximum acceptable sampling can be found from geometrical con- 
siderations. Since the signals of interest are confined to angles near the nor- 
mal to the scan plane, the sample spacing should be such that neighboring 
array elements are 180° out of phase when the incident angle approaches 
the boundary of the region of interest. This is the first null of the array 
pattern. The sample spacing is given by: 

A 

ASm ““ 2 cos(| + 0) 

where 6 is the incidence angle. The resulting grating lobes are located at 
incidence angles 4>\ 


, tt , n\ 

<f> = g ± arccos( — ); 


n = 1,2,3, . 


For the PMTC specifications, the target zone width is 20 feet at a 
distance of 59 feet from the scan plane aperture. Thus, the maximum 
angle of interest is roughly 10 degrees from the scan plane normal and the 
maximum acceptable sample spacing is 2.8A. 
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Table 1: Position Accuracy Requirements 


f [GHz] 

A [m] 

A/10 [mm] 

A/360 [mm] 

0.5 


60.00 

0.1667 

1.0 


30.00 

0.0833 

18.0 


1.666 

0.0463 

35.0 


0.857 

0.0238 

phase error 

360° 

36° 

1° 


2.4 Position Accuracy 

Near field position must be known to within a small fraction of a wave- 
length to obtain accurate field measurements. For system frequencies of 1 
to 18 GHz, the minimum wavelength is 1.667 cm. For a maximum 1° error 
in the out of scan-plane direction, the z' coordinate must have 0.046 mm 
accuracy. For a maximum 10% A error and the additional requirement of 
1.667 mm accuracy across the diagonals, the allowable error in x' and y' 
directions can be found to be 1.6/\/2 = 1.17 mm. This assumes indepen- 
dence of ( x',y ') and z' errors. These results are summarized in Table 1. 
Note that the highest frequency of interest is the limiting factor. 

Note that phase measurement errors and cable bending errors can be 
modeled as z' coordinate errors. Bearden and Dugenske [5] report that 
Georgia Tech has a 1° accuracy over their scan plane. Yaghjian [8] notes 
that position error may be reduced by scanning along both vertical and 
horizontal lines and appropriately averaging the two data sets. This how- 
ever, doubles data set measurement times. An important concern is re- 
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positioning accuracy, necessary for background subtraction algorithms (if 
needed) and system reliability. All errors must be contained such that the 
total effective error is within the tolerances stated above. 

2.5 Probe Characteristics 

Johnson et.al. [1 pp.1671] suggest the following probe characteristics: (1) Small 
aperture for point measurement and broad beam width. (2) Polarization 
should be accurate and stable. (3) Probes the size of A/2 dipoles or smaller 
are needed for accurate measurements but probes with more directivity 
than A/2 dipoles produce large measurement errors. 

To reduce the time required for each data set scan, different scanning 
antennas can be used. Cown and Ryan in [9] suggest the use of arrays of 
modulated scattering elements or scanning antenna arrays to accomplish 
one or two of the dimensions of the scan plane. The data measurement times 
for these antennas are shown to be reduced by a factor of approximately 
100 in comparison to scanning with a single probe antenna. [9 pp.580 Fig. 8] 

2.6 Scanner to Test Target Separation Dis- 
tance 

2.6.1 Evanescent waves 

An evanescent wave’s constant phase plane is not parallel to its direction 
of arrival. Thus their propagation constant can be greater than k Q , the 
free-space propagation constant along a given direction [6]. Incident an- 
gles ignored in the sample spacing computations may also correspond to 
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evanescent mode propagation constants. Due to these abnormalities, the 
evanescent wave effects should be minimized. One solution is to place the 
antenna array far enough away from the scatterer such that the evanescent 
waves are attenuated by the distance. Joy and Paris [3] formulated the 
minimum attenuation of the evanescent waves at a distance, NA. 

a min = 8.7N\(kl - kl)i [dB] 

where 

ki < < k\ + 

k\ = fc 2 sin 2 9 cos 2 ^ 
fc 2 = fc 2 sin 2 6 sin 2 <f> 

Note that sin 2 9 is greater than 1. 

2.6.2 Multiple Interactions 

Multiple interactions are all the possible reflection paths from the target 
to the scan plane via walls, support structures, etcetera. Multiple inter- 
actions can be reduced by increasing the scan surface to target distance, 
using efficient absorber, decreasing probe size and averaging the results of 
scans fractions of a wavelength apart in the z' direction. To experimentally 
estimate the effect of multiple bounce interactions, Yaghjian [8] suggests 
varying the spacing between the scan plane and target as a function of the 
wavelength. Amplitude variations of period one half wavelength will be 
caused primarily by double bounce reflections. 
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2.7 Conclusions 


We have mentioned some of the important parameters to consider in the 
design of a near field scanner. Joy [4] makes an exhaustive list of param- 
eters, some of the more notable being: (1) Record the position, ampli- 
tude and phase of the probe to implement error compensation algorithms. 
(2) Smooth motion of the probe antenna will minimize vibration. (3) Mo- 
tion can induce Doppler effects on the frequencies. (4) Allow for measure- 
ment of both polarizations. (5) Test facility temperature should be held 
constant with minimal gradients. (6) Teflon has a phase anomaly at room 
temperature. 
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Chapter 3 

Planar Scanner Simulation 

3.1 Introduction 

We modeled the planar near-field scanner as an antenna array in the trans- 
mitting mode, generating a plane wave region in the target zone. By reci- 
procity and the duality of a scanning probe and an antenna array, the planar 
near-field scanner receive pattern is similar to the transmission pattern of 
the antenna array. Thus, use of the transmitting viewpoint to design the 
planar scanner is justified. 

3.2 Theory 

The simulation is accomplished by a point by point application of Maxwell’s 
wave equation solution in free space for spherical wave propagation. Each 
field point in the quiet zone is the resultant sum of contributions from each 
source element antenna in the scan plane. The equation used to model 
the antenna array system is described in Figure 3. Normalization of the 
E n sources in the scan plane by N, the number of array elements, gives a 


13 



constant output power independent of the numbers of array elements. 

The amplitude and phase taper functions were chosen to simulate the 
edge treatment effects of a compact range reflector with rolled edges. In 
a rolled edge compact range, the reflector geometry transitions from a 
parabolic surface, which produces a constant phase and amplitude reflected 
signal, to a rolled edge, which produces an increasing phase delay and a 
decreasing amplitude. In the planar scanner taper function a similar tran- 
sition from constant phase and amplitude to phase and amplitude tapering 
occurs. 

This taper function, T(n), is used to spatially filter the measured data 
to suppress measurement sidelobes. The taper function is illustrated in 
Figure 4, The amplitude taper implemented is a cosine squared rolloff 
which characteristically has low side lobes, and the phase taper is a cosine 
function drop off with a frequency independent phase delay characteristic. 
A given taper function is denoted as the percentage of the scan plane width, 
x% of W and the maximum phase delay (MPD) in degrees, for example 
14%, 90°. 

3.3 Tests Implemented 

The figures to be discussed are one-dimensional cuts through what is, in 
practice, a two-dimensional pattern. These patterns are the E-field am- 
plitude and phase that result in the quiet zone from the excitation of the 
array antenna in the scan plane. The normalized E-field amplitude is given 
in units of dB/m which means 201og 10 (V7m) where V/m is the measured 
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Err, 

— 

N E pjkRnm 

E " T(n) P(<? nm ) 

n— 1 Jy 

where 

T{n) 

= 

Taper function 

P(0 nm) 

= 

Scan plane element pattern function 

L 

= 

angle of line joining 11 th and m th elements 

k 

= 

free space propagation constant 

Rnm 

= 

distance between n lh and m th elements 

N 

= 

Number of elements in scan plane 

D 

= 

distance between scan plane and quiet zone 

AS 

= 

Scan plane element spacing 


Figure 3: Simulation Equation Development 




PHASE OF T (n) (DEG) MAGNITUDE OF T (n) (V/m) 



SCAN PLANE POSITION (n) 



Figure 4: Visualization of Taper Function 
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E-field amplitude. The E-field phase is referenced to the scan plane and is 
given in degrees. The constant parameters of the simulations presented are: 

(1) The scan plane array element pattern function is set to 1 for all 0 nm , 
sinndating an isotropic radiator. (2) The scan plane width, W is nominally 
8.5 m or approximately 40 ft. The physical dimensions used were closely 
matched to the design specifications of the PMTC Milicon-P986 Bistatic 
RCS measurement range. See Figure 5. 

We investigated the effects of (1) amplitude and phase taper variations, 

(2) increasing the element spacing, A S', beyond A/2, (3) changing the dis- 
tance between the test zone and scan plane, and (4) the response of the 
model to frequencies outside the design frequency band. In each case, the 
aperture size attained, W qz , is discussed. Comparison is made with the 
characteristics of the compact range reflector with rolled edges at the Ohio 
State University as reported by Burnside et.al. [11]. 

3.4 Results 

3.4.1 Taper Effects 

The first parameter investigated was the amount of amplitude taper which 
gave the largest quiet zone for a 1 dB, 10° approximation to a plane wave. 
The frequency was chosen to be 1 GHz as a worst case where the ripple 
would be greatest. The phase taper was held constant at 90° while the 
amplitude taper was varied. The results are shown in Figure 6 and tabu- 
lated in Table 2. Note from Table 2 that the phase ripple has a minimum 
whereas the amplitude ripple stays constant. 
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Figure 5: Bistatic RCS Measurement Facility with Planar Scanner 
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Table 2: Scan Plane Array Antenna Amplitude Taper Response, 90 


Taper [% 

Ripple [dB] 

Ripple [deg] 

W qz [m] 

Size Ratio 

12 

1.2 

15 

5.6 

1.52 

13 

1.1 

13 

5.4 

1.57 

14 

1.0 

11 

5.2 

1.63 

16 

1.0 

12 

5.2 

1.63 


Table 3: Scan Plane Array Antenna Phase Taper Response, 12% 


Taper [deg] 

Ripple [dB] 

Ripple [deg] 

W QZ [m] 

Size Ratio 

45 

1.5 

12 

5.0 

1.70 

90 

1.2 

15 

5.6 

1.52 

135 

1.1 

13 

5.4 

1.57 

180 

1.0 

12 

5.4 

1.57 


Next the amplitude taper was held constant at 12%, and the phase taper 
was varied to try to bring the ripple to within the 1 dB, 10° specifications. 
The results of this are shown in Figure 7 and tabulated in Table 3. From 
Tables 2 and 3, note that similar responses were obtained for tapers of 
12%, 180° and 14%, 90°. We chose the 14%, 90° taper for the tests to 
follow, leaving more room for variation of the phase taper if necessary. 

3.4.2 Sample Spacing 

Figures 8, 9 and 10 illustrate the effect of element spacing greater than 
A/2 at 1, 10 and 18 GHz respectively. The amplitude and phase ripple are 
summarized in Table 4. It is apparent from these figures that a sample 
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Normalized E-Field Phase (deg.) Normalized E-Field (dB/m) 








Table 4: Scan Plane Array Antenna Sample Spacing Variation Response 


f [GHz] 

Spacing [A] 

Ripple [dB] 

Ripple [deg] 

W qt [m] 

Size Ratio 

1.0 

0.5 1 

1 

11 

5.2 

1.63 

1.0 

1 

12 

5.6 

1.52 

1.5 

2 


5.6 

1.52 

2.0 

5 1 

24 

4.0 

2.13 

2.5 

8 

67 

2.2 

3.86 

10.0 

0.5 

1 

3 

6.4 

1.33 

1.0 

1 

3 

6.4 

1.33 

1.5 

1 

3 

6.4 

1.33 

1.75 

1 

3 

6.4 

1.33 

2.0 

1 

9 

5.8 

1.47 

2.5 

1 

12 

2.6 

3.27 

18.0 

r 0.5 

1 

3 

6.6 

1.29 

1.0 

1 

3 

6.6 

1.29 

1.5 

1 

3 

6.6 

1.29 

1.75 

1 

3 

6.6 

1.29 

2.0 

1 

9 

6.2 

1.37 

2.5 

1 

12 

2.8 

3.04 


spacing of 1.75A is satisfactory for the higher frequencies, but at 1 GHz, 
1A is the largest satisfactory spacing. It can be calculated from Table 4 
that an element spacing of 2.6 cm in the scan plane is acceptable for all 
frequencies below 18 GHz. The reason that 2A spacing has a poor quiet 
zone width is the close proximity of the grating lobes to the main beam of 
the antenna array in the scan plane. 


22 



















































Normalized E-Field Phase (deg.) Normalized E-Field (dB/m) 


1GHz D-12.49m Taper 14%,90°AS=0.5 1.0 1.5a 



-5 -4 -3 -2 -1 0 


-- 0.5A 

- 1.0A 

- 1.5A 

3 4 


Quiet Zone Position (m) 


2 3 

Quiet Zone Position (m) 


Figure 8: Normalized E-Field Amplitude and Phase, Comparison of Differ 
ent Sample Spacings at 1 GHz 



Normalized E — Field Phase (deg.) Normalized E Field (dB/m) 







Normalized E-Field Phase (deg.) Normalized E-Field (dB/m) 


18GHz D=12.49m Taper 14%,90°AS=0.5 1.75 2 A 


0.5A 

1.7 5 A 

2 A 


-3-2-10 1 2 3 

Quiet Zone Position (m) 


0.5A 

1.75A 

2A 


- 3 - 2-10 1 2 

Quiet Zone Position (m) 


Figure 10: Normalized E-Field Amplitude and Phase, Comparison of Dif 
ferent Sample Spacings at 18 GHz 







3.4.3 Aperture Size 

For a worst case of 1.0A spacing at 1 GHz, our design gives a scan-plane 
to target-zone size ratio of 1.5, see Table 4. This is an improvement over 
the scanners mentioned in the literature which have size ratios only as low 
as 1.7, and it is much smaller than the OSU compact range reflector ratio. 
For the OSU compact range reflector at 1 GHz, the amplitude ripple is 
1 dB and the phase ripple is 10° over a 2.7 m quiet zone height [11 pi 79 
Fig. 10(a)]. Dividing the nominal height of the compact range reflector (5.8 
m) by the quiet zone height results in a scan-plane to target-zone ratio of 
2.1 . 


3.4.4 Down Range Variation 

Figures 11 and 12 illustrate the performance over the front, middle and back 
of the quiet zone for frequencies of 1 and 18 GHz respectively. At 1 GHz, 
even though the responses are all within a 1 dB npple, the distance change 
causes differing ripple functions as the range increases. This indicates a z 
dependence of the plane wave approximation which must be minimized. For 
18 GHz, the ripple functions are all well suppressed and no z' dependence 
is observed. The results of these test are given in Table 5. 

3.4.5 Outside Design Frequency Band 

The last two figures, 13 and 14, show the response of the scanner array 
to the frequencies 0.5 and 35 GHz which are outside the design frequency 
band of 1 to 18 GHz. The ripple is tabulated in Table 3. Note that a 
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Normalized E-Field Phase (deg.) Normalized E-Field (dB/m) 


1GHz D=10.37 12.49 14.61m Taper 14%,90°AS=0.5X 




Quiet Zone Position (m) 


Figure 11: Normalized E- Field Amplitude and Phase, 1 GHz Down Range 
Response Comparison 
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Normalized E-Field Phase (deg.) Normalized E-Field (dB/m) 











Table 5: Scan Plane Array Antenna Down Range Performance 


f [GHz] 

Range [m] 

Ripple [dB] 

Ripple [deg] 

W qz [m] 

Size Ratio 

1.0 

10.37 

1 

10 

6.0 

1.42 

12.49 

1 

11 

5.2 

1.63 

14.61 

1 

14 

5.4 

1.57 

18.0 

10.37 

1 

3 

6.6 

1.29 

12.49 

1 

3 

6.6 

1.29 

14.61 

1 

3 

6.6 

1.29 


Table 6: Scan Plane Array Antenna Performance Outside Design Band 


f [GHz] 

Spacing [A] 

In— toUEM 


w„ 

m 


Size Ratio 

0.5 


14,90 

1,12 

3.8 

2.24 

35.0 

■M 

22,135 

1,3 

6.6 

1.29 


minimum scan-plane to target-zone ratio of 2.2 at 0.5 GHz still compares 
well with the OSU compact range ratio of 2.1 at 1 GHz. 
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0.5GHz D=12.49m Taper 22%,135°AS=0.5 A 


- 3 - 2-10 1 2 3 

Quiet Zone Position (m) 


Quiet Zone Position (m) 


Figure 13: Normalized E-Field Amplitude and Phase, 0.5 GHz Respon 
Below Design Frequency Band 








35GHz D=12.49m Taper 14%,90°AS=1.75A 
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Figure 14: Normalized E-Field Amplitude and Phase, 35 GHz Response: 
Above Design Frequency Band 






Chapter 4 
Conclusions 


The one dimensional results given above are directly applicable to the two 
dimensional case with possibly slight modifications to the windowing func- 
tions. The sample spacing tests suggest the use of a constant separation 
of the elements in the array (or the probe sample spacing), making the 
sample spacing independent of the frequency. The maximum theoretical 
spacing of 2.8A however, is not attainable due to the close proximity of the 
grating lobes. The variable distance tests indicate a need to quantify the 
z' dependence effect and reduce it for low frequencies. The tests outside 
the design band indicate a possibility to use these frequencies with some 
modification to the filtering or tapering functions and a possible decrease in 
target-zone width. An advantage of this planar scanning system is that per- 
formance modifications can be accomplished without hardware adjustment. 
Modifications or upgrades are done only on the processing algorithms. 
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Chapter 5 
Future Work 

5.1 Scanner Construction 

To complete this project, we will build a prototype scanner to test the 
theoretical computations presented and to expand to the two dimensional 
case. We will test the necessity of optical /laser position tracking systems 
using computer error modeling as well as direct manipulation of the scanner. 
We will try to quantify the z* dependence of the fields at low frequencies 
and will try to find or develop an optimal filtering function of the scan 
plane antenna elements. We will extend the scope of the scanning process 
to theoretically evaluate the behavior of a scanner based on a scanning 
parabolic strip. 

5.2 RCS Calibration Technique 

Once the scanner has been designed and is operational, it can be used to 
measure calibrated bistatic RCS. The process of measuring bistatic RCS is 
similar to that for monostatic RCS. In contrast to a compact range reflector 
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system, planar scanner data sets need spatial filtering before being used to 
find the bistatic RCS. The measurement steps are as follows: (1) Scan a 
calibration target, usually a sphere. (2) Scan the background. (3) Vector 
subtract the background from the calibration target data set. (4) Apply 
theoretical bistatic RCS of sphere to normalize the experimental value ob- 
tained. 
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Appendix A 
C Program Listing 
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(include <stdio.h> 

(include <math.h> 

(define MAX_IO_PTS 2000 
(define MAX 10 COL 4 

(include <ll5ERT: (REUNING.C. INCLUDE ]DEF .H> 


/* This program computes the fields in the 

/A**************************************** 

/* Write data to FileName, w/header 

/A**************************************** 

void WriteFile{ DataBlk , numpts , numcol , title 
double DataBlk [ MAX_IO_COL ][ MAX_I0 
int numpts, 

numcol ; 
char *title, 

♦axistitles , 

♦FileName ; 

{ int LCV , LCV2 ; 

FILE *outf ileptr ; 


quiet zone from an aperture */ 

******************♦*♦♦♦******♦**♦*♦*/ 

V 

*****************♦♦*♦*******♦****♦*♦/ 
ax i sti ties, FileName ) 

/* data to print */ 
/* array access */ 
/♦ number of columns */ 
/* data graph titles */ 
/* axis title array */ 
/* output file name */ 

/♦ loop control variables */ 
/* output file pointer */ 


if ( ( outf ileptr - fopen( rileName , "v" ) ) !- NULL) /* file open error? */ 

* printf ( "\nWriting %d columns to ts .\n\n" .numcol .rileName ) ? /* info */ 


printf < "(s", title ) ; /* print the titles */ 

fprintf(outfileptr,"%s", title)? 

printf { "%s" , axistitles) ? /* print the axis titles */ 

fprintf(outfileptr,"%s" , axistitles ) ; 


for (LCV - 0? LCV <« numpts-1; LCV++ ) 


/* print the data */ 


for ( LCV 2 - 0; LCV2 o numcol-1? LCV2++) 

fprintf {outf ileptr, "Ilf ", DataBlk [ LCV2 )[ LCV]) ; 
f printf (outf ileptr ,"\n" ) ? 


) 

fclose(outfileptr) ? 

printf ( "Write Completed\n " ) i 


/* close output file */ 
/* user info output ♦/ 


} 

else 

printf ( "Could not open filename - 


/♦ error handeler */ 

% s .\n* , FileName ) ; 


} /* end procedure writeFile */ 


/* Main ***************************♦♦***♦♦♦♦****************/ 

main( ) 


/♦ declare variables */ 


char outf ilename[60] , 

title 1 120 ] , 
mainti tie [ 50 ] , 
subtitle! 70 ] , 
axistitles! 120 ] , 
xaxistitlel 35 ] , 
ylaxistitlel 35 ] , 
y2axistitle[ 35 j ; 

double freq, 

lambda , 


/♦ filename variable */ 
/♦ data graph titles */ 
/♦ main title string */ 
/* subtitle string V 
/* axis title array V 
/* x title both graphs ♦/ 
/* yl axis title V 
/* y2 axis title V 


/* operating frequency in GHz */ 
/* wavelength in cm */ 
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int 


double 


double 


double 

double 


complex 


di st , 
apln, 

exhalfapln, 
apsp, 
qzln , 
qzsp, 
rolloff , 
maxphsdlay ; 


/* separation of aperture and quiet zone in cm */ 

/* aperture length in cm */ 
/* half* the aperture lenath in cm */ 
/* aperture spacing in cm */ 
/* quiet zone length in cm */ 
/* quiet zone spacing in cm */ 
/* roll-off in amplitude/phase taper in cm */ 
/* max phase taper delay in degrees */ 


numap , 
numqz , 
testid , 
lev, lcv2 ; 


/* number points in aperture */ 
/* number points in quiet zone */ 
/* test id number */ 
/* loop control variables */ 


apem[MAX_IO_PTS] , 
apept MAX_IO_PTS ] , 
apex [ MAX_IO_PTS ] , 
apey [ KAX_IO_PTS j ; 


/* aperture E field magnitude in V/m */ 
/* aperture E field phase in radians */ 
/* aperture E field X coordinate in cm */ 
/* aperture E field Y coordinate in cm */ 


qzem[MAX_IO_PTS] , 
qzep[MAX_IO_PTS) , 
qzex[MAX_IO_PTS] , 
qzey [ MAX_IO_PTS ] ? 


/* quiet zone E field magnitude in v/m */ 
/* quiet zone E field phase in radians V 
/* quiet zone E field X coordinate in cm */ 
/* quiet zone E field Y" coordinate in cm */ 


block [ MAX_IO_COL H MAX_IO_PTS ] ; /* data block for file I/O*/ 


r /* length between current point antennas in cm */ 
T h e ta /* angle of exit from aperture in radians */ 
PAT . ' /* element pattern: function of Theta */ 


Eap[HAX_IO_PTS] , 
Eq Z ( MAX 10 PTS ) , 
Etemp, Curlle; 


/* aperture E field */ 
/* quiet zone E field */ 
/* temporary variables to use in summation */ 


/* function declarations *******/ 


complex C add(),C sub{ ) , 

C_rect { ) ,£_ini t ( ) , 
C_mult ( ) ,C_exp( ) ? 


/* complex addition and subtraction */ 
/* convert to rect format and initialization */ 
/* complex multiplication and exp function */ 


double get double(), 
powT ) , 
atan2 ( ) , 
cos( ) , 
logl 0 { ) , 
C_mag( ) , 

C_ang ( ) ? 


/* function to get user input */ 
/* exponentiation function */ 
/* arctangent function */ 
/* cosine function */ 

/* log base 10 function */ 

/* returns the magnitude of a complex number */ 

/* returns the arguement of a complex number */ 


int get_int { ) ; 


/* user input function */ 


/********************♦******************/ 
/* input the general setup of the^test^V 


testid 
freq - 


. qet int ( "\nEnter the test ID number — >",0,0); 
get_double ( "\nEnter the FREQUENCY in GHz -->”,0.5,40.0)) 


lambda ■ 30.0/freq; 


/* in cm */ 


di st 
apsp 

apln 


jet double( "\nEnter the SEPARATION in ■ — 0.0,40.0*100.0; 
jet^double ( "\nAperture ELEHENT SPACING in WAVELENGTHS — > , 
).0,5.0)*lambdaj . . . „ 

jet double ( "\nHaximum Aperture WIDTH in ■ •-> ,0.0, 20.0) *100.0, 


numap « ( apln/apsp ) +1 ; 
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apln - ( numap- 1 )*apsp; /* the exact length */ 

exhalfapln - apln/2.0 + apsp; 

i f ( numap >- MAX_I0_PTS-1 ) pr intf ( "\nNumber of element* exceed* array size. ); 

qzln - get double { "\nMaximum Quiet Zone WIDTH in « — >" , 0 . 0 , 20 . 0 ) *100 . 0 > 
qzsp - qzln/100.0? /* in cm. ensure* 100 points in the quiet zone */ 


numqz - { qzln/qzsp ) +1 ; 
qzln - ( numqz-1 ) *qzsp? 
i f ( numqz >- MAX 10 FTS-1 ) 


/* the exact length */ 
printf ( "\nNumber of elements exceeds array size."); 


rolloff - (get double ( "\nROLLOFF percentage — > " , 0 . 0 , 50 . 0 )/100 . 0 ) *apln ; 
maxphsdlay - get_double ( "\nMax . PHASE LAG in degrees — > " , 0 . 0 , 360 . 0 ) ? 
printf ( "Enter a main title. \n"); 

gets ( mainti tie ) ; /* discard the return from the previous input */ 

gets(maintitle) ; 


/* set up the position of aperture and quiet zone arrays^*/ 

f or ( lcv-0 ; lcv<-numap-l ; I cv-m- ) 

{ 

apex( lev ] -0 . 0 ; 

apey( lev ]- ( apln/2 . 0 )-{ (lev)* apsp) ; 

} 


f or ( lcv2-0 ; lcv2<-numqz-l ; lcv2++ ) 

{ 

qzex(lcv2J-dist; 

qzey[lcv2)-{ qzln/2 . 0 ) - ( (lcv2)*qzsp) ? 

} 


/**************************************** * ** */ 

/* set up the values of the aperture arrays */ 

/+ a*********#********************************/ 

f or( lcv-0 ; lcv< -numap- 1 ; lcv++ ) 

{ /* Amplitude taper in V/a. can be any function of position (lev) */ 

if (apey(lcv) <- exhalfapln && 

apey(lcv) >- exhalfapln - rolloff it 
rolloff I - 0.0) 

apemt lev)- ( pow( cos (( PI/{ 2* rolloff ))*( apey [ lev J-exhalf apln+rolloff ) ) 
, 2 . 0 ) /numap) ; 

else if(apey[lcv] >- -exhalfapln £& 

apeyilevj <- -exhalfapln + rolloff &£ 
rolloff l- 0.0) 


{ 

apem[ lev 1- 
) 

else 

{ 

apem[ lev ]- 
) 


( pow ( cos ( (Pl/(2*rolloff) )*(apey[ lev J+exhalf apln) + 
(PI/2 . 0 ) ) , 2 . 0 ) /numap ) ? 


( 1 .0/numap) ; 


/* Phase taper in radians can be any function of position */ 
i f ( apey [ lev ] <- exhalfapln tfc apeyllcv) >- exhalfapln - rolloff) 

apep[ lev)- maxphsdlay* ( PI/180 . 0 )*( sin ( ( PI/( 2 . 0* rolloff ) )* 

( apey [ lev J-exhalf apln+rolloff )+Pl/2 * 0 )-l ) ; 

else ^ if { apey [ lev ) >- -exhalfapln && apey(lcv) <- -exhalfapln + rolloff) 

apept lev ) - maxphsdlay* ( PI/180 . 0 )*( sin( (PI/(2.0*roHoff ) ) * 

(apey[ lev j+exhalf apln ) )-l ) j 
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else 


) 


{ 

apepllcv]- 0.0; 

} 

/* combine into a complex number */ 
Eapt lcv]-C_rect (C_init(apem[ lev] ,apep! lev ]* 180. 0/3. 14 1592653 59 ) ) ; 


printf ( "\nProcessing . 

/* calculate the quiet zone fields */ 

for(lcv2«0 ; lcv2<-numqz-l ; lcv2++ ) 

' Etemp-C_init(0. 0,0.0) ; 

f or ( lcv-0 ; lcv<-numap-l ? lcv++ ) 


/* user information */ 


{ 


R - pow( ( pow( dist , 2 . 0 )+pow( apev [ lev ]-qzey[ lcv2 ] , 2 . 0 ) ) , 0 . 5 ) ; 

Theta - atan2{ (apey[lcv]-qzey[lcv2) ) ,pow( diet, 2.0) ) ; 

PAT - l; /* any function of Theta «/ 

CurEle - C mult(C mult(Eap[ lev] , C init( (PAT/F) ,0.0) ), 

C~exp( C_init ( 0 . 0 ,-7*3 .14159265359/lambda*R) ) 


); 


) 


Etemp - C_add( Etemp, CurEle ) ; 

1 

Eqz] lcv2 ] - Etemp; 

qzem[ lcv2 ] - 20 . 0*logl0 ( C_mag( Eqz [ lcv2 ] ) ) ; 
qzep( lcv2 j - C_ang( Eqz [ lcv2 ] ) ; 


/* in dB/m */ 
/* in degrees */ 


/* output y coordinate, magnitude t phase^*/ 


"QZ%03d.dat",testid) ; 


/* plot title generation */ 


sprint f (outfilename, ' 
sprintf ( subtitle , 

" %gGHz D-%gm Taper %g \% ,%g S-%g \\", 

f req,di st/100. 0,rollof f/apln*100 . 0 , maxphsdlay, apsp/lambda) ; 
sprintf (title, "%s |td\n%s\n" .rnainti tie , testid, subtitle ) ; 
sprintf (xaxistitle,"[l]Quiet Zone Position (m)")j 
sprintf(ylaxistitle,"[2 ] Normalized E-Field (dB/m)"); 
sprintf(y2axistitle,"[3 ] Normalized E-Field Phase <d*g-> )» 

sprintf (axisti ties, "ts\n%s\n%s\n", xaxistitle, ylaxistitle,y2axisti tie) ; 


printf ( "\nStoring. . .\n" ) ; 

for ( lcv-0 ;lcv<-numqz-l ; lcv++) /* put data into the data block */ 

block] 0] l lcv]-qzey[ lev ]/l 00.0; 
block] 1 j ] lev j -qzem] lev ] ; 
block] 2] ] lcv]-qzep] lev] ; 

Wr iteFile( block, numqz, 3, title, axisti ties, outf ilename ) ; 
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